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GLOBAL EXISTENCE OF WEAK SOLUTIONS TO THE 
THREE-DIMENSIONAL PRANDTL EQUATIONS WITH A 
SPECIAL STRUCTURE 

C.-J. LIU, Y.-G. WANG, AND T. YANG 


Abstract. The global existence of weak solutions to the three space dimen¬ 
sional Prandtl equations is studied under some constraint on its structure. 
This is a continuation of our recent study on the local existence of classical 
solutions with the same structure condition. It reveals the sufficiency of the 
monotonicity condition on one component of the tangential velocity field and 
the favorable condition on pressure in the same direction that leads to global 
existence of weak solutions. This generalizes the result obtained by Xin-Zhang 
on the two-dimensional Prandtl equations to the three-dimensional setting. 


Contents 


1. Introduction 

2. The proof of the main result 

2.1. Porous medium-type equation 

2.2. Transport equation 

2.3. Proof of Theorem ll.il 
References 


1 

4 

5 

_8 

n 

ig 


1. Introduction 

Consider the initial boundary value problem for the Prandtl boundary layer 
equations in three space variables, 


( 1 . 1 ) 


'dtu + {udx + vdy + wdz)u -I- dxP = d^u, in Q, 

dtv + (udx + vdy + wdz)v + dyp = d^v, in Q, 

< dxU -\- dyV + dzW = 0, in Q, 

= 0, lim {u,v) = (U{t,x,y),V{t,x,y)), 

z—¥-\-(yc> 

^(u,'u)|t =0 = {uoix,y,z),vo{x,y,z)), 


where Q = {t > 0,(x,y) € D,z > 0} with a fixed D C (U(t,x,y),V(t,x,y)) 
and p{t, X, y) are the traces of the tangential velocity field and the pressure of the 
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Euler flow on the boundary {2: = 0}. Note that the traces satisfy 


dtU + UdxU + VdylJ + dxP = 0, 
dtV + UdxV + VdyV + dyp = 0. 


Despite of its importance in physics, there are very few mathematical results 
on the Prandtl equations in three space variables. In fact, except the recent work 
[5] about the classical solution with special structure and those in the analytic 
framework [Ilia, most of the mathemtical studies on this foundamental system 
in boundary layer theory are limited to the problem in two space dimensions, cf. 
[niiiiTiEiiioiiisiiii] and the references therein. 

Recently in we obtain the local well-posedness of classical solutions to the 
problem (HU under some constraint on the structure of the solution, in order to 
avoid the appearance of secondary flow (i) in boundary layers. Precisely, assuming 
that for the Euler flow given in (II.2p . U{t, x, y) > 0, in the class of boundary layers 
that the direction of tangential velocity field is invariant in the normal variable z, 
and the x—component of velocity uit, x, y, z) is strictly increasing in z, dzU > 0, in 
0 we have constructed a classical solution to the problem (HH]), and it is linearly 
stable with respect to any three dimensional perturbation. In the class of this 
special structure, the solution of (HU takes the form: 


(1.3) {u{t, X, y, z),k{x, y)u{t, x, y, z),w{t, x, y, z )), 

where the function k(x,y) satisfies the following condition (H): 

(HI) the function k depending only on {x,y) satisfies the inviscid Burgers equa¬ 
tion in D, 


(1.4) kx + kky = 0; 

(H2) the outer Euler flow 

(U(t, x, y), k{x, y)U(t, x, y), 0,p{t, x, y)) 
with U{t,x,y) > 0, satisfies that from the system (11.21) . 

f dtU + UdxU + kUdyU + dxP = 0, 
l^dyp - kdxP = 0 . 

Moreover, the authors recently observed in [B] that for the shear flow z),v‘‘{t, z), 0) 
of the three-dimensional Prandtl equations, the special solution structure HU is 
the only stable case. 

Under the above assumption HU of special solution structure, the original prob¬ 
lem HU of three dimensional Prandtl equations is reduced to the following one for 
two unknown functions (u,w)\ 


(1.5) 


where dQ- 


dtu + (udx + kudy -I- wdz)u — d^u = —dxP, in Q, 
dxU + dy(ku) -\- dzW = 0, in Q, 

< 

mU=o = wU=o= 0, lim u = U(t,x,y), 

2 —)- + 00 

MdQ- = Wi(i,a:,?/, z), u|t=o = uo{x,y,z), 

= (0,oo) X 7_ X R+, with 

7- = {ix,y) G dD\ {l,k{x,y)) ■n{x,y) <0}, 
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and n{x, y) is the unit outward normal vector of D at (x, y) S dD. 
For this reduced problem under the assumption that 


( 1 . 6 ) 


dzUQ > 0, dzUi > 0, for z >0, 


in the class of dzU > 0, we apply the method developed by Oleinik m for two 
dimensional Prandtl equations. Precisely, by the Crocco transformation, 


^ = ri = y, C 


u{t,x,y,z) 

U(t,x,y) 




dzu{t,x,y,z) 

U{t,x,y) 


the problem dLSl) becomes the following initial boundary value problem, 

(1.7) 

{ L{W) = dtW + CU{d^ + kdrj)W + AdcW + BW-W'^d^W = Q, in fir, 
1F|C=1 = 0, lF9clF|c=o = t- 
W\r_=W,it,^,yX) = ^, 

W\t=o = WoX,V,0 = ^, 


where 

^T = {{tX,vX)\o<t<T, X,v)eD, 0<C<1}, 

T- = {(t,^,?7,c)l 0 < t < T, X,T]) e 7 _, 0 < C < 1}, 

and 

(1.8) ^ = -C(l - C)^ - (1 - C")^, B=^+(:{U,+ kUy) - dyk ■ Ct/. 

For the problem (13 of the degenerate parabolic equation, we established local 
existence of classical solutions in [ 5 . 

As a continuation of the paper [^, the purpose of this paper is to prove the 
global in time existence of a weak solution to the problem uni for data satisfying 
dm) and the favourable pressure condition: 


(1.9) Pi{tX:V) < 0, fort>0, X,f])eD. 

For this, we will adopt the approach introduced by Xin-Zhang in [1^ for the the 
two-dimensional Prandtl equations to the three-dimensional setting. As observed 
in m, the main motivation of introducing the favourate condition on pressure is 
to avoid the separation of boundary layers. 

For completeness, the definition of weak solutions to (II3 is given as follows. 


Definition 1.1. A function W{tX,V,0 G L°°(^0,T\BV{V,)'^ for some T > 0 is 

called a weak solution of the problem ini) int < T, if the following conditions hold: 
(i) There exists a positive constant C such that 


c-\i-0 < w{tx,vX) < 0(1-0, y{tx,y,0&^T. 


(ii) W satisfies the first equation and the initial boundary conditions of uni) in 
the weak sense: 

pT 


/o jvl 


+ CiUkip)y + {Aip)(; - Bfj - Wip(^(^'^df,d'qdC,dt 


In VFo 


■'ip\t=odf,dT]dC - 


/o JD 


U ' W 


l(-oXdr/dt 


/ 7 . 

/o ay_ Jo 


Wi 


kndt^dldt, 
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for any test function S C°°{flT) satisfying 

■0 = 0, at i = T or (^, rj) € 7 +; 0^ = 0 at C = 0. 

Here, 

1+ = {{,x,y)&dD\ {l,k{x,y))-n{x,y)>Q}, 

and the function kn = ( 1 , k) ■ ft with ft being the unit outward normal vector of 

The main result of this paper can be stated as follows. 

Theorem 1.1. For the problem ini), and any T > 0, assume that k G C^(D), 
U € C^((0,T) X D), px G C^((0,T) X D) satisfy (11.91) . and the initial boundary 
data Wo G G C^(r_) satisfy 

(1.10) Co\l-0<Wo,Wi<Co{l-0, 

for a positive constant Cq- Then, there exists a weak solution W{t,f^,r],f) G 
L°^^0,T;HV(fl)^ to the problem (11.71) in the sense of Definition 1.1. 


2. The proof of the main result 

Following the approach introduced in El: a viscous splitting method is used 
to construct a sequence of approximate solutions to the problem dni). Precisely, 
divide the time interval [0,T] into n equal sub-intervals: 

T 

0 = to < < • • ■ < tn-i <tn = T, ti+i — ti = —, for 0 < t < n — 1. 

n 

First, in the time step [0, ti] or [F, ti+i] for an even i, we construct the approximate 
solution by solving the following initial boundary value problem for a porous media- 
type equation: 

{ ^dtW — W'^d'^W + Ad^W bW = 0, in (0,ti] x Tt or (ti,ti+i] x 12, 
lF|t=o = Wo, or W\t=u = W{ti,f,r],C) (given in the previous step), 
iF|c=i = o, Wclc=o = t; 

and in the time step [ti,ti+i\ for an odd i, we construct the approximate solution 
by solving the following problem for a transport equation: 

{ \dtW + CU{d^ + kdr,)W + {B- b)W = 0, in (0,0+i] x 12, 

W\t=ti = W{ti,^,r],f) (given in the previous step), 
lF|r_ = IFi. 

Here, the coefhcient function b in will be chosen later to satisfy the boundary 
condition VF|r_ = Wi for the solution of (12.1|) . What is needed is to prove that 
the function W constructed in the time interval [0,r] is uniformly bounded in n 
and has a uniform total variation with respect to the spatial variables rj and f. 
This implies that as n —> 00 , the limit function of the approximate solutions W 
constructed in (I2T]1-(I^ is a weak solution to the problem dnj. The proof is 
divided into the following subsections. 


PRANDTL EQUATIONS IN THREE SPACE VARIABLES 


5 


2.1. Porous medium-type equation. In this subsection, consider the following 
problem for a porous medium-type equation 

{ dtW - W'^dlW + AdcW + bW = 0, for 0 < t < T, (^, ??, C) e 
W\t=o = Wo{^,V,0, 

< 0, Vb|c=i = 0. 


In order to match the boundary condition W\r_ = Wi, by observing that Wi > 0 

for 0 < C < 1 and 


Wi = 0(1 - C), as C ^ 1, 


we set 

(2.4) + Ajr^ ■ 

where /(^, r/) is a non-negative smooth function defined on the closure of the domain 
D satisfying = 1. By the formulations in (11.81) : 


A = -(1-C) 


^^ + (1 + ^)^ 
^ u ^ u 


and the assumption given in Theorem 11.11 there exists a positive constant Mq 
depending on the parameters of |L3, such that for the function b(t,^,r],(^) given 
in (1^ . 


(2-5) \\b(t,^,riX)\\w^'°°{nT) ^ ^0- 

The problem (12.31) can be viewed as a one space dimensional problem by regarding 
variables ^ and rj as parameters. 

Note that the equation in (12.31) is degenerate on the boundary = 1}. As [2], 
consider the following uniformly approximated parabolic problem: 

r dtW, - (w^ + €)d^^w, + Ad^w, + bw, = 0, for t > 0, (C, V.0 e 

( 2 . 6 ) \w,\t^o = Wo>Q, 

[lT,dclT.|c=o = ^<0, TT,|c=i = 0, 

for a positive constant e > 0. It is known that problem p.6|l has a unique smooth 
solution. After getting some uniform bounds of We with respect to e, we can obtain 
a solution to the problem (12.31) by taking e —>■ 0 in (12.6p . In fact, we have 


Theorem 2.1. Under the assumption of Theorem 11.11 the problem (lOl) has a 
unique solution W € BV(0,T;il) for any fixed T > 0. Moreover, W has the 
following properties: 

(1) there exists a positive constant jd, depending on ||1To||l°o and the L°°-norm of 
the parameters of (lOl) . such that 

(2.7) ^oe-^'V < W{t, e, 7?, C) < Cie^^^l - C), 

where 


(/j = 62 ? sin|(l - C), 6»o = min{Wo/v?}, 
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(2) there exists a positive constant C 2 , depending on ||T4^||loo, ||9(jWo||l“ OLnd the 
-norm of the parameters of such that 

(2.9) |W^cl<C^2; 


(3) 

for any t > 0, 



(2.10) 

[\wc{t,f,pX)\dC< 1 

JO j{ 

' \Woxitv, 

C)\dC + W{t,f,r],0)-Wo{^,V,0) 

and 




(2.11) 


Jo ^ 

-C)'dC + C3t(i + ^Vo.cMc) 


^0 JO 

W^s) 


Here, the positive constant C 3 depends on ||VF||j;,oo and the C^-norm of the pa¬ 
rameters of (ESI). Also, Wri and Wt satisfy similar estimates as (12.1111 by simply 
replacing the partial derivative with respect to f in (12.1111 by the partial derivatives 
with respect to rj and t, respectively. 

Proof. To prove the inequality (12.7L we first show W\(^=q > 0 by using the fa¬ 
vorable condition on the pressure function. For this, we can assume that ^ < 0 
because we can replace it by ^ — i5 for some constant d > 0 and then let <5 —>■ 0. 
Then, we want to show that 

(2.12) lFe|c=o >0, Ve > 0 

holds for the problem (12.611 . Otherwise, by using the continuity of and lFo|(^=o > 
0, there exist eo > 0 and a point P on {C = 0}, such that Weo|p = 0. That is, 
attains its minimum at P, which implies that d(A/Vf^\p > 0. But from the 
boundary condition of (1^ on C = 0, we have 

lc=o = = 2^ < 0, 

which is a contradiction. Hence, we obtain (I2.12I1 . which implies that VF|c=o > 0 
for the problem (12.311 by letting e —>■ 0. 

Now, combining lF|(^=o > 0 with the boundary condition WWq\q=q = ^ < 0, 
we obtain < 0- The rest of proof for (12.711 is similar to that of Lemma 4.2 

in [14] so that we omit the details. Moreover, there is a constant m > 0, such that 

(2.13) lF|^=o > m, We\c=o > m, Ve > 0 
hold. 

(2) We now turn to the estimate (|2.9p . Consider the problem (12.611 for W^, and 
the corresponding problem for i9^1Fe as follows, 

- 2W,5<;H4a^W, -h (A^ -h b)d^W, = -b^W,,, 
dcW,lt=o=Wox, 

W,5cVF.|c=o = ^, a2lF,|c=i=0, 

where the operator 

L^ = dt- {W^ + €)-dl + A- dc- 

Here, note that 9|We|^=i = 0 because We|^=i = 0 and = 0 in the first 

equation in (12.611 . 











PRANDTL EQUATIONS IN THREE SPACE VARIABLES 


7 


Set V = dcWe — a( with a being a constant to be chosen later. It follows that 

L'^(V) - 2W^{V + aC)V^ + iA^+b- 2aWe)V 

= i2a^C - -aA- aCiA^ +b)^Y. 

From the first step, Y is bounded provided that a is bounded. Let Vi = e~^*V 
with /3 > 0 being sufficiently large satisfying 


P + A(; + b-2aW^ > 1, or /3 > \\At; + b - 2aW^\\L-- + 1. 

From (12.1511 and (I2.14L we have 

(2.16) L^{Vi) - 2W,{V + aC)(Vi)c + {(3 + A(; + b - 2aW,)Vi = 
with the following initial and boundary conditions 

(2.17) foi|t=o = lF,foi|c=o = e-''*^, dcVi\c=i =-ae-^K 

Firstly, note that for an arbitrarily fixed constant a > 0, from (I2.17|) . we have 
i9(;Fi|^=i < 0. Also, from (I2.13P and the relation in (12.171) 

VF,Vi|c=o = e-''*^<0, 

it implies that Vi does not attain its positive maximum on ^ = 1 or ^ = 0. Then, if 
Vi attains its positive maximum in the interior of or when t = T, we have from 
(12.161) that 


Vi < maxje ^*Y} < ||y||j;,oo. 

If Vi achieves its positive maximum when t = 0, it follows that 
Vi < max|lFo,c - aC} < llbFo.clU”- 

Therefore, we conclude that Vi < max|||F||L“, ||||lo°|, which implies that 

(2.18) dcWe < a + e^*max|||y||L~, ||WoxI|l“|- 

Secondly, for an arbitrarily fixed constant a < 0, by considering the possible 
negative minimal points of Vi on Ht, similar to the above argument, we have 

foi>-max{||r|U~,|lbFo.clU~,llf7^|^^J|L->} 


which implies that 


(2.19) 


diYVe >a- e^‘max|||y||Loo, |11 Fo,c||loo, 




Hence, combining (12.181) with (12.1911 and letting a —)> 0 yield that 


\dtA¥^\ < max III y Hi-, ||lFo,clk“> 


Px 1 

u w. 



where /3 > ||A,^ + 6||i- + 1 and Y = —b^W^- Thus, we obtain (12.91) as e —>■ 0. 

(3) The proofs of (I2.10|l and (12.111) are similar to those given in Lemmas 4.6 
and 4.7 of respectively. And the proof for the uniqueness of solution to the 
problem (12.311 is similar to that of Theorem 4.1 in El- Thus, we omit the detail 
for brevity and this completes the proof of the theorem. □ 
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2.2. Transport equation. In this section, we will study the problem of transport 
equation ( 12 . 2 |) for t S with i being odd, that is, we consider the problem, 

r \dtW + C,U{d(^ + kdr,)W + hiW = 0 , in x fl, 

( 2 . 20 ) 

[lT|r_ =Wi, 
where the function 

(2.21) bi(t, C V, C) = {B- b){t, r], () 

with functions B and b being given in (II. 8 p and (12.411 . respectively. 

For any fixed (t, ? 7 , C) G (^i, ^i+i] x 11, the characteristics of the equation (j2.20l) 

passing through this point is denoted by: 

(s, 72(s;t,C,^,C), c), s€{ti,U+i] 

with 7 i and 72 being determined by 

= 2 C-c/(s, 7 i(s;^,C, 0,72(5; 

72 ( 5 ; 1,0 7,0 = 2C-(fcti)( 5 , 71 ( 5 ; 1,0 7, 0 , 72 ( 5 ; 1,0 7 , 0 ), 
(7i,72)(1;1,07,0 = (07). 

For simplicity of notations, in the following we will also use abbreviations 71 ( 5 ) 
and 72 ( 5 ) to represent 7i(5;i,^,?7,0 and 72 ( 5 ;!,07,0 respectively, when without 
confusion. 

Combining (12.2211 with the property + kkrj = 0, we have 

^fc(7i(5),72(5)) = 0, 

which implies that 

(2.23) A:(7i(5),72(5)) = A:(C,7), Vs S (t,,ti+i]. 

Then, from (12.2211 and (12.2311 it follows 

(2.24) 72 ( 5 ) = fc(07)-7i(5) + 7-fc(07)0 
and 7 i(s) is given by 

^2 25 ) / Vi(®) = K-U (s, 7i(s), k{i, rf) • 71 ( 5 ) + 7 - 1^(0 7)C) , 

l7i(l) = 0 

Observe that 7 ^( 5 ) > 0, and the projection of this characteristic on the (^, ry)-plane 
is a straight line passing through (^, 7 ) with slope /c(^, 7 ). Moreover, the function 
A:(07) remains constant along this line. 

Note that the solution of the above problem (12.2511 exists and is unique when 
the function is Lipschitz in (5, 7 ). The solution W of the problem (j2.20ll 

is represented by 

(2.26) = W^(s, 7 i(s), 72 (s),c) exp|-y 61 (s, 71 ( 5 ), 72 ( 5 ),c)ds|. 

Denote by 

(2.27) 

1 *( 1 ,?, 7 ,C) = inf |l e [U,t] : Vs S {i,t], (71(5; 1 , 7 , C), 72(5; 1 , 7, O) S nj. 
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Note that if t* > U, then 

e dD. 

Denote by 

Qi = e {ti,u+i] X n 

and 

Q 2 = [(t,^,r],C) € {ti,t,+i]x n : > U, (7i(i*).72(i*)) e 7-}, 

where 7 _ is the closure of 7 _ = {(^,??) G dD : {l,k{^,ri)) ■ < 0} on the 

boundary dD. 

To study the estimate of the solution to problem (12.2011 , we first give the following 
proposition for the representation of the solution. 


Proposition 2.2. For the problem (12.201) . we have 

(tij X D — U Q2: 

and at any (t, C) ^ (ti,ti+i] x fl, the solution of (12.201) can be represented by 

(2.28) 

W{t,^,r]X) 


fT(^t,_i, 7 i(ti_i), 72 (t*-i),c) -expl - 5 i(s, 7 i(s), 72 (s),C)ds|, in Qi, 
W^i(^t*,7i(t*),72(t*),c) - expl - 6i(s,7i(s),72(s),C)ds|, in Qa- 


Proof. It suffices to show that {ti,ti+i] x fl = Qi U Q 2 , because the formulation 
(|2.28|1 follows immediately from (12.261) . We divide the proof into the following two 
parts. 

(1) Firstly, we claim that 


{t, 77 , 0) G Qi, for all t G (t,, ti+i], (C, ry) G D. 

Indeed, from (12.241) and (12.251) it follows that 7 ((s) = 0, and ( 7 i(s), 72 ( 5 )) = (C, 17 ), 
which implies that t*{t, rj, 0) = ti by using the definition (12.271) . 

(2) Next, we will prove that if 

^)\Qu C > 0 , 

then {t,f,r],C.) ^ Q 2 - Indeed, for such point {t,^,riX), there is t*{t,^,r],() > ti 
such that 

{C,v*) = {ii{n,72ir)) &dD. 

We then need to show that P = (C*, i?*) G 7 _. 

From (12.241) and (|2.25|1 . we have 

(2.29) 7 ((s)> 0 , 7 i(s) > r. Vs G (t*,t], 

and 

(2.30) ^^(^s) = k{C,V*){li{s)-C)+V*, VsG(r,t]. 

Without loss of generality, assume that in a neighborhood Ps of the point P = 
(^*, 77 *) in the {f, ? 7 )-plane, the boundary of D is represented by a smooth function 
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77 = /(C), SO that T] > /(C) when (C, r]) G Dr\ Ps- Then, the outward normal vector 
at the point ^C, /(C)) is given by 

(2.31) n(C, /(C)) = , ^ (/'(C), - 1 ). 

\/l + (/'(C)) 

By the definition of t*, we have there exists a eo > 0 such that 
72 ( 5 ) >/(7i(s)), Vse(r,r + eo]. 

From (12.291) and (I2.30L we know that there exists a constant ei > 0 such that 

(2.32) vCe (C*,C*+ei]. 

If P does not belong to 7 _, then there exists a di < d such that 

(1, *(C, ??)) • ^(C, 7) > 0 , V(C, 77) e 511 n Ps^ 

which implies from (12.311) that there exists a £2 > 0, satisfying 

( 1 , fc(C, /(O)) • (/'(C), -1) >0, VC G (C* - £ 2 ,C* + £ 2 ). 

This is, 


(2.33) /'(C)>fc(C,/(C)), VCG (C*-e2,C*+e2). 


For a fixed Co G (C*, C* + £3) with £3 < min{£i, £2}, consider the function 

^(c) = /(c)-/(Co)-fc(c,/(c))(c-Co), CG r,co]. 

Since + kk^ =0 from dm, it follows that 


i"'(C) = /'(C)-fc(C,/(C))- 
= /'(C)-fc(C,/(C))- 
= /'(C)-fc(e,/(e))' 


%(C,/(C))+/'(C)fc,(C,/(C))] -(C-Co) 

-fc(C,/(C))+/'(C)]A,(C,/(C))-(C-Co) 

i-fc,(C,/(C))-(C-Co)l, 


and then, from (12.331) and the fact that is bounded, we have 
(2.34) i"'(C)>0, VCG(C*,Co) 


provided that £3 is sufficiently small. Therefore, (12.341) implies that 


0 = F(Co) > F{C) =V*- /(Co) - HC,ri*){C - Co), 

which is a contradiction to (12.321) by letting C = Co in (12.321) . Hence, we have 
P G 7 _, which implies that (t,C,??, C) G <32, and this completes the proof of the 
proposition. 

□ 
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2.3. Proof of Theorem II.IL Based on the results obtained in the above two 
subsections, we will give the proof of Theorem 1 1.1 1 in this subsection. Before it, the 
following lemmas and propositions are needed. 


Lemma 2.3. Let W be the approximate solution constructed by and (1^ . 

Then 

(2.35) |tT(t,e,i7,C)l <5ie"^‘(l-C), 

where 

Cl = maxICi, = max|Mi, \\B - 6 ||l=o| 

with positive constants Ci and Mi being given in (12.81) . Moreover, there exists /3 
depending only on ||Wo||loo, |11L"i||c2 and the parameters of problem (11.71) . such that 


(2.36) 


W{t,f,,r],C) > 9oe ^*(p, 


where g) is given in 


and 9 q = min 


Wo 
V ' 


Wi ] 

V /■ 


Proof. When 0 < t < ti, the estimates (|2.35|) and (12.361) follow from (12.7p in 
Theorem 12.11 immediately. Assume that (12.3511 holds for 0 < t < ti with i > 1, and 
consider the case for tt <t < t^+i. If i is even, from (12.7L we have 

\W{t)\ < maxjll^lU^, - C) < Cie^^\l - C), 

by using the induction hypothesis. If i is odd, from (I2.28P and (12.211) it follows 
\W{t)\ < \W{ti)\ ■ < Cie^i*(l - 0, 


by using the induction hypothesis again, or 

\W{t)\ < Wi ■ < Cie^i‘(l - 0- 

Thus, we conclude the estimate (12.351) . 

Next, suppose that (12.361) holds for 0 < t < with i > 1. Then if i is odd, from 
(IQ) in Theorem 12.ll we have that there exists ^ depending on Ci, ||&||l°° and the 
parameters in the problem (E3 such that 

W > min{^^}e-«‘-‘^V > ^oe-^V- 

If i is odd, the estimate (12.361) is a direct consequence of the expression (|2.26p in 
Proposition [221 provided that /3 > ||i? — Thus, we complete the proof. □ 


Remark 2.4. From Lemma \2.3\ and by virtue of (1^ . we find that there exists a 
constant Co, depending only on ||Wo||loo, HtTilIca and the parameters in the problem 
(ini), and satisfying Co > Co with the constant Co being given in (11.101) . such that 

(2.37) Co-i(l - C) < W{t,^,p,0 < Co(l - C). 

Now, we study the estimate of the first order derivatives of the approximate 
solution with respect to the spatial variables for obtaining the uniform estimate on 
the total variation of the solution. Before it, we give the following two propositions 
for the problem (12. 2p of transport equation. 
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Proposition 2,5. For the problem (I2.2p . there exists a constant C 4 depending on 
the domain D, the eonstant C'o given in Ij2.37ll . ||W^i||ci o,nd the estimates of 
the parameter in the problem (j2.2p . sueh that for all t G [U^ti+i] and f G (0,1), 


(2.38) 


ID 




W^it) 

+ C4(t — ti) C 4 


ti JD 


|W^g(5)| + |W^,(5)| 

W^{s) 


(1 — O^d^dpds. 


Proof. For the problem (12.2|) . we know that satisfies 

(2.39) 

+ caj [c'( j(=)£] + + awu^), - (b - 

Taking p.37[l into account, we multiply the above equation ()2.39|1 by {l — ()‘^signW^ 
or (1 — C)^ ■ , and integrate the resulting equation over D with respect to {f, rj), 

to obtain that 

+ Cd, (1 - Cfd^dp - C(fcC/),H(l - Cfd^dp 

< \m)^\\L^ • /^ ^(1 - cyd^dp + \\B- 611^0.. H(i - cYd^dp 

+JJiB-b),^l^\d^dp. 


(2.40) 


From (12.2F we obtain that on the boundary 7 _, 

CC/(A:,5, + kndn)W\^_ = -dtWi - {B - b)\^_ ■ VFi, 

which implies 

(2.41) CUkr^dnW\^_ = -dtWi -{B- b%^ ■ IFi - CU\^_ ■ krO^Wi ^ 62 . 

Obviously, there exist two bounded functions ai{f,p) and a2{Cv) defined on the 
boundary 7 _ such that 

= aidn + a 2 dr, on 7 _. 

Thus, from (j2.41|l one has 

(2.42) CUkr,W^\^_ = 0162 + a2CU\^_ ■ kndrWi ^ 63 - 
Hence, it follows that by virtue of (12.421) . 

(2.43) 


.|lFd 


iw^d, 


f D 


dfdp= I CU-^{1- 0 il,k)-ndl 


> 


I CUkJ-^ii-C?di = - I&3I 


dD 

(1-C)" 


Wf 


-di 


>-(Oo)dl& 3 ||L~ 0 ( 7 -), 

where l{p-) is the length of 7 _ and the positive constant C'o is given in (ll.lOF 
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By using (j2.37p . it follows 

(2.44) JjiB- 6)5< CoKB - 6)d|L~ • 3 ( 0 ), 

where S{D) is the area of the domain D. 

Plugging (12.431) and (I2.44|) into (12.401) , we obtain that there exists a constant C 4 
depending on D, Cq, ||fPi||ci and the estimates of the parameter in the problem 
(12.21) . such that 


Tt 


'D 


W^{t) 


(1 - Crdidr] <C4 + C4 




ID 


W^{t) 


(1 - Cfd^dr]. 


Similarly, we can obtain another constant C 4 such that 




D 


W^{t) 


(1 - C,fd£,d'q <C4 + C4 


f \Wdt)\ + \W,{t)\ 

Id W^t) 


{l-Cfd^dr]. 


By letting €4 = 04 + 04 , we have that from (12.451) and (12.461) . 
(2.47) 


± [ mit)\ + \wm 

dt Jd W^t) 


(1 - Qfd^dri <04 + 04 [ 

Jd 


\wdt)\ + \w^it)\ 

W^{t) 


{l-Cfd^dr]. 


Then, integrating the above inequality (I2.47P over {ti , t) gives the estimate ()2.38l) 
immediately, and we complete the proof of the proposition. □ 


Proposition 2.6. For the problem (12.2|) . there exists a constant O 5 depending on 
the domain D, the constant Co given in (j2.37l) . ||VPi||ci and the O^ estimates of 
the parameter in the problem ()2.2p . such that for t G [Ujti+i] and C. £ (0,1), 

(2.48) 

J \Wi;{t)\dfdri < J \W(;{U)\d^dr] + 05 {t-U) 

Proof. From the problem (12.21) . we know that Wc^ satisfies 

(2.49) dtWc, + CU{d^ + kd^)Wc + U{W^ + kW^) + {B - b)W(; = -{B - b)c^W. 
Multiplying the above equation (12.491) by signWQ or > and integrating over D 
with respect to (^,??), it follows that 

(2.50) 
d f 


dt 


'D 


mmig + C(9«(C/|lFcl) + dr,{kU\Wc\))d^dp - C 


ID 


\W^\d£,dT] + ||fcl7||j 


< l|t/||L= 

+ \\(B-b)cWh^-SiD). 

As in (I2.43|) . we have 


\Wr,\d^dr] 


\B-> 


ID 


U^ + {kU)r, ■ \Wi;\dfdr] 
\W(;\d^dr] 


ID 


(2.51) 


Id 


• C|lTcl • (l,fc) d^dr, > -||fc„C/VFi,cl|Loo(^_) • Ih-)- 
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Obviously, from the bounded estimate (12.371) for W one has 


(2.52) 

/ mdfdr,<{Cof 

JD 


and 



(2.53) 

[ \Wr,\d^drj < {Cof 

JD 


Plugging (12.511) . (12.521) and (12.531) into (12.501) . we obtain that there exists a constant 
C 5 depending on D,Co, ||fPi||ci and the estimates of the parameter in problem 
(EH), such that 


(2.54) ^ J^lWc(tMdv<C5 + C5 [llPcWI 




W^(t) 




d^drj. 


Then, integrating the above inequality (12.541) over the estimate (12.481) in the 

proposition follows immediately. □ 


Remark 2.7. Similar to the above proposition, one ean show that Wt satisfies an 
estimate similar to (j2.48l) in Proposition IS.fi] by replacing the partial derivative in 
C by that in t. 

It is ready to give the estimate of the first order derivatives of the approximate 
solution W constructed in (I2.1l) - (j2.2p with respect to the spatial variables. 


Lemma 2.8. For any fixed T > 0, let W{t, (0 < t < T) be the approximate 

solution of dnp constructed in Then, there exists a constant M > 0, 

depending on T and the domain D, the constant Cq given in (j2.37l) . ||Wo|1l=o, 
|| 1 Ti||c 3 and the estimates of the parameter in the problem (II.7p . such that for 
all t G [0, T], we have 


(2.55) 


/n 


\Wc\ + m\ + \W^\ it,-)d^dr,dC 


< m(i + • J^i\Wox\ + |W"o,d + \WoJ)d^dpd(:) . 


Proof. The proof is divided into the following three steps. 

(1) When t G {ti,ti+i] for even i,W is determined by the initial boundary value 
problem p.ll) for a porous medium-type equation. From Theorem 12.11 and the 
boundedness (j2.37p of W, we obtain that for t G (ti,ti+i], 

(2.56) 


lWt;l(t,-)d^dr/d( < / lWi;l(p,-)d^dr/d( 


W{t,f,,r],0) -W{P,£,,t], 0) d£,drj. 


ID ■- 


Moreover, there exists a constant C 3 depending on the domain D, the constant Cq 
given in (I2.37L ||lFi||c3 and the estimates of the parameter in problem (II.7L 
such that 

-m\ + \w^\ 


In 


< 


1F2 

m\ + \Wr, 


(1 - C) (t. ■)df,dr]dC 


In 
f C 3 


n 


(1 - O'] {U,-)d^dr,dC + ^ 3(1 + ^ \Wc\{U, fid^dgdc) ■ {t - U) 
(1 — O' (S; fid^drjdCds. 


1 F2 




M/2 
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which implies that by using the Gronwall inequality, 

- 0^1 {tr)d^dr]dC 


/n 






(2.57) 






(1-C) {U,-)d^dridC 


< • / 

Jn 

+ (gCalt-q) _ 1) . I'l + ^ \Wc\iU,-)d^drjdc) . 


Combining (I2.56P with ()2.57|) . it follows that 

fibbcl + - 0"] (t, ■)d^drjdC 

' Vcl + - C)"l iu, ■)d^dr,dC 

d^dr]. 


(2.58) 


Jn 

< gC 3 (t-ti) . 


Jn 

+ e^^d-u) 


ID 




C=o 


(2) When t S {ti, ti+i] for odd i, we obtain W by the problem (12.21) for a transport 
equation. From Propositions 12.5112.61 and the estimate (12.371) , it follows that there 
exists a constant Cq, depending on the domain D, the constant Cq given in (I2.37L 
llbFilIca and the estimates of the parameter in problem (II.7L such that for 
t € {ti , tj+l], 


In 


< 


|lFcl + 

|Wc| 


m + \w^\ 


(1-C)^ {t,-)d^dr]dC 


in 




+ Ce 


ti J Q, ^ 


|lFc| + 


(1-C) {ti,-)d^dridC+ Ce{t-ti) 
(1 — C)^ {s, ■{d^drjdCds, 






which implies that by using the Gronwall inequality, 


(2.59) 


|Wc| 

in 






(1-C)" {t,-)d^dr,dC 


in 


l^cl + - 0"] {U,-)d^drjd(: + - 1. 


(3) On the other hand, when i is odd we have that from Proposition [5^ 


(2.60) 


W(ti+i,C,??,0) = W{ti,^,r),0) • exp| - ^ {B - 6)(s, C, 0)ds|. 


By combining (12.581) , (12.591) and (12.601) . and letting C 7 = maxICa, Ce}, we 
obtain that for any t G {ti,ti+i], 


(2.61) 


|VFcl + 

In L 

< gC 7 (t-ti) . 




1 F2 

IVFcl 


(1 - C) (i, ■)d^dr]dC 

(1 -C)"l {U,-)d^dr]dC 




Jn 

+ gC7(t-q) 
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where 


Id 


C=o 


d^drj 


I W(t,-)■ exp (S - 6)(s,-)rfs|j 


for even i, and Gi{t) = 0 for odd i. Hence, 

(2.62) G,(f) < 2 ||TH||loo • S{D), Vf S [0,T]. 
Therefore, (I2.6ip implies that 

(2.63) 


in 


llTcl 






(1-C) {t,-)d^dridC 


< ■ f [|lTo,cl + - C)"] d^drjdC + - 1 + G,(t) + F,(t), 


where 


F,(t)=^e^^(‘-‘^+^)G,(t,+i) 


3=0 


) gC^TCi-tj+i) 


0<j<i—2, j'.odd 

Note that Wo|^=o > 0, then 

(2.64) 


f W{tj+i,-)■ exp^- {B - b){s,-)ds'^ 


C=o 


d^drj. 


m) 


< e 


2Ct^ 


fT|c=oL • [w"o(-) • exp { - J\b - 6 )(s, ■)ds} 


C=o 


d^drj 


41-1 


f r / T ri2k+i V 

^ lT(f 2 fe,-)- (l-exp{- 2 G 7 -/ (H - 6)(s, ^ds} j 

k—l ^ ^ d t2k 


C=o 


d^drj 


< 


lT|c=o ■ S{D) ■ + I), 


where [|] is the largest integer less than or equal to |, and 


1-1 


^=E 


:,C’r{t — t2k) 


k^l 


1 — exp{— 


rt2k+l r- 


'^2k 


2Cr + {B- b){s,^,r],0) ds} 


By choosing a constant G 7 satisfying that 2 G 7 > 
141-1 


{B-b)\c^o 


L°°{D) 

we obtain 




,C7^(i-2fc) . 


('t2k + l r 


(2.65) 


f41-l 


'i2k 


2 G 7 + (H - 6)(s,^, 77,0) ds||L=(D) 


< V . 4c-^_ < 

n 
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Plugging p.62|l . (12.6411 and (12.6511 into (I2.63L it follows that there exists a positive 
constant Cg, depending on T,S{D), ||iP||Loo and Cy, such that 


JD 
< e 


liPcI 

Cjt f 






(1-C) it,-)d^dr] 


ID >- 




K‘ 


d^dr] + Cg, 


from which the estimate (j2.55p follows immediately by using (j2.37p . Thus, we 
complete the proof of the lemma. □ 

We are now ready to give the proof of the existence of weak solution as follows. 

Proof of Theorem 11.11 Denote by ly, C) the approximate solution con¬ 
structed in (IOl-(l^. Let Y be the dual space of First, we claim that 

'(1-C)^' 


( 2 . 66 ) 


\ J 


L2(o,T;y) 


< Cg 


with a constant Cg independent of n. Indeed, if i is even, we have that from 
/■‘i+i /ii _yi2. 2 r 1 5 1 

= L.'r Ja K'"’” + - C)‘ 


dt 


fti+i p 


< Cio 




with a uniform constant Cio depending only on D, ||4l|l(;yi and ||6||icx,. 

If i is odd, it follows that by using (1^ . 

/•*i + l /D— 2 /■*i + l r 1 R _ /in 


dt 


<Ci 




1^^-^ 




with a uniform constant Cu depending only in D, ||?7||ci and \\B — &||loo. Thus, 
by using (I2.37L the estimate (12.6611 follows immediately. 

Next, from (12.3711 and Lemma [2.81 we know that 

(1-C)^ 




< Cl 2 , 


L“(0,T;lvm(n)) 

where C 12 is a positive constant independent of n. Hence, by using the Lions- 
Aubin Lemma (see [2] for instance), we conclude that | | is compact in 

L^((0,T) X D). Therefore, we may assume that 

(i-C)^ (1-C)^ 




IT 


-, in D{{0,T)xn), 


and then 


IT" 


IT, a.e. in (0, T) x D. 
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In particular, 

in L^(0,T) X fi). 

Thus, for any ip £ (7“ ((0, T) x fl), we have that 
(2.67) 


i'.odd 


fn L 


- - A( - )c + CUi - )e + CUki -)„ + 






2b-B 

IT" 


ipd^dridC,ds. 


From Theorem o and Remark \2l[ and by integrating by parts in the above 
equality (12.671) it yields that. 


( 2 . 68 ) 


>0 Jn 


IT" 


i>t + CiUip)^ + C{Uk'ip)n - {B - b)'ip 


d^drjdCds 


= 2E 


i'.odd *■ 


£ + aukid\ + (26 


- B 


d^drjdC,ds. 


Letting n —>■ oo in p.68|l and noting that ti — ti-i = ^, we have 


pT 


'ifjt + C(^^)c + Q{Ukil))n — {B — b)ilj d^d7]dC,ds 

+ ^ ( “ + C{Uip)i + CiUk'4’)n + {2b-B 


d^dridC,ds, 


Jo Jn 
which implies 
(2.69) 

I + CiUkid), + (Ai6)c - - Wi,cc}d^drjd(:ds = 0. 

Therefore, from ()2.69p we know that the function IT satisfies the equation of prob¬ 
lem (EZl) in the sense of distribution. Moreover, we can obtain that IT satisfies the 
estimate (I2.55P by letting n —>■ oo, which implies that 

IT G L^{0,T;BV{n)). 

It remains to verify that IT satisfies the initial and boundary conditions given 
in (II.7|) . It follows from Lemma E751 that. 


M,W(t,i,n.O = 0. a.e.m(0 ,T)xa 

We can verify the other boundary conditions in dni) for IT in the sense of distri¬ 
bution, respectively, through similar process as above to show that IT satisfies the 
equation of dni in the sense of distribution. For example, we have the following 
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equality holds 


/o Jn 


1 

Wl 


i’t + + C{Ukip)rf — {B — b)il) d^dr]d(ds 




+ {2b-B 


d^drid(ds — 


1 


IdU’ W 


\c=od^drids, 


for any 

,peC^{{ 0 ,T)xDx{-l,l)) 

with = 0. Therefore, we complete the proof of Theorem I l.li 


□ 
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